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ABSTRACT 

Context. The extragalactic distance scale builds directly on the Cepheid Period-Luminosity (PL) relation as delineated by the sample 
of Cepheids in the Large Magellanic Cloud (LMC). However, the LMC is a dwarf irregular galaxy, quite different from the massive 
spiral galaxies used for calibrating the extragalactic distance scale. Recent investigations suggest that not only the zero-point but also 
the slope of the Milky Way PL relation differ significantly from that of the LMC, casting doubts on the universality of the Cepheid 
PL relation. 

Aims. We want to make a differential comparison of the PL relations in the two galaxies by delineating the PL relations using the 
same method, the infrared surface brightness method (IRSB), and the same precepts. 

Methods. The IRSB method is a Baade-Wesselink type method to determine individual distances to Cepheids. We apply a newly 
revised calibration of the method as described in an accompanying paper (Paper I) to 36 LMC and five SMC Cepheids and delineate 
new PL relations in the V, I,J,&K bands as well as in the Wesenheit indices in the optical and near-IR. 

Results. We present 509 new and accurate radial velocity measurements for a sample of 22 LMC Cepheids, enlarging our earlier 
sample of 14 stars to include 36 LMC Cepheids. The new calibration of the IRSB method is directly tied to the recent HST parallax 
measurements to ten Milky Way Cepheids, and we find a LMC barycenter distance modulus of 1 8.45 ± 0.04 (random error only) from 
the 36 individual LMC Cepheid distances. In the J, K bands we find identical slopes for the LMC and Milky Way PL relations and 
only a weak metallicity effect on the zero points (consistent with a zero effect), metal poor stars being fainter. In the optical we find 
the Milky Way slopes are slightly shallower than the LMC slopes (but again consistent with no difference in the slopes) and small 
effects on the zero points. However, the important Wesenheit index in V, (V - 1) shows a metallicity effect on the slope and on the zero 
point which is likely to be significant. 

Conclusions. We find a significant metallicity effect on the W vl index y(W V i ) = -0.23 ± 0.10 mag dex~'as well as an effect on the 
slope. The A"-band PL relation on the other hand is found to be an excellent extragalactic standard candle being metallicity insensitive 
in both slope and zero-point and at the same time being reddening insensitive and showing the least internal dispersion. 

Key words. Stars: variables: Cepheids, Stars: fundamental parameters, Stars: distances, Magellanic Clouds, distance scale 



■ 1 . Introduction of the quest for the Hubble constant. A dissenting view can be 

found e.g. in Sandage et al. (2009) and references therein. 

The Cepheid Period-Luminosity (PL-) Relation is fundamental 

to the calibration of the extra-galactic distance scale and thus to ^ value of the PL relatlon rests wlth lts universality, in 
the determination of the Hubble constant. Modern reviews on Particular that the PL relation slope and zero points are inde- 
the calibration of the Cepheid distance scale can be found in e.g. P endent of metallicity. The zero point has long been suggested 
Freedman and Madore (1991), Fouque et al. (2003), Sandage t0 be metallicity dependent and the HST Key Project on the 
and Tammann (200^, Fouque et al. £2007), and Barnes G3E> Extra g alactlc Dlstance Scale (Freedman et al. [2001) corrected 
while Freedman and Madore <2SD review the present status for thls based on the empirical studies available atfliat time (e.g. 

Freedman and Madore 11 9901 Sasselov et al. 1997, Kochanek 
119971 Kennicutt et al. 1998 ) with an estimated effect in the 



* Based on observations collected at the European Organisation W w index of ~ 02 ± 02 ma § dex ~ in the sense that metal-poor 

for Astronomical Research in the Southern Hemisphere, Chile, Cepheids are fainter than metal-rich Cepheids. However, the size 

Programme-IDs 076-C.0158, 078.D-0299, & 080.D-0318 and ev en the sign of the effect is still disputed. Udalski et al. 

** Table 2 is only available in electronic form at the CDS (120011 ) found no effect on the slope nor on the zero-point when 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via comparing the metal -poor Cepheids of IC16 13 with those in the 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ Magellanic Clouds based on TRGB distances to those galax- 
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ies. On the other hand Romaniello et al. (|2008) found a signif- 
icant effect with the opposite sign when comparing Milky Way 
and Magellanic Cloud Cepheids with individual spectroscopic 
metallicity determinations. Sakai et al. (120041 1 and Storm et al. 
(2004), found an effect of similar size and sign as that adopted 
by the HST key project, but the latter study also showed a sig- 
nificant difference in the slopes of the Milky Way and LMC re- 
lations, as also found by Sandage, Tammann, and Reindl (2004 ). 
Recently Bono et al. (2010) have argued for no significant effect 
in the Wvi and Wjk index on the basis of empirical as well as 
theoretical arguments. 

Observationally, we are in the dilemma of either using a large 
sample of Cepheids, like in the Large Magellanic Cloud (LMC), 
which constrains very well the slope of the relation for a low 
metallicity sample of stars and which leaves the zero point to 
be determined from secondary indicators, or of using direct geo- 
metric distances (parallaxes) to a handful of nearby, solar metal- 
licity, Milky Way Cepheids which constrain well the zero point 
of the PL relation but which do not constrain the slope very well. 

Baade-Wesselink type methods which use the pulsational 
properties of the Cepheids to determine direct distances to in- 
dividual Cepheids promise to resolve this dilemma by yield- 
ing direct individual distances to a large sample of Milky Way 
and Magellanic Cloud Cepheids spanning a significant range of 
metallicities. In particular, the near-infrared surface brightness 
(IRSB) method (Fouque & Gieren, 119971 ). a near-infrared vari- 
ant of the Barnes-Evans method (Barnes and Evans, [19761 ) shows 
great promise in achieving this goal as it is insensitive to red- 
dening errors. The technique has also been shown to be inde- 
pendent of the metallicity of the Cepheids (Storm et al. 2004). 
A few years ago the IRSB method was re-calibrated using in- 
terferometrically measured, phase-resolved angular diameters of 
Cepheids by Kervella et al. (2004|. This was an extremely im- 
portant achievement as it proved that the pulsating Cepheid vari- 
ables indeed obey the same surface brightness-colour relation as 
the stable yellow giant stars which were used to construct the 
surface brightness-colour relation adopted in the original cali- 
bration of the technique (Fouque & Gieren, 1997 ). The improved 
version of the IRSB technique was then applied to Milky Way 
Cepheids by Storm et al. (120041 ) and Groenewegen (2008), and 
for the first time to extragalactic Cepheids by Storm et al. (2004 ) 
(SMC) and Gieren et al. (I2005al ) (LMC). 

In our earlier application of the method to Milky Way 
Cepheids (Storm et al. 2004), we found that the slopes of the 
resulting PL relations in all optical and near-infrared bands were 
significantly steeper than the those directly observed in the Large 
Magellanic Cloud by the OGLE Project and by Persson et al. 
(2004), thus challenging the universality of the PL relation. 
Gieren et al. (2005a) then analyzed thirteen LMC Cepheids, for 
which the data required for the IRSB analysis were available at 
the time, in an identical fashion and found that the PL relation 
slopes in the LMC were very similar to the one obtained for 
the Milky Way sample, but different from the observed apparent 
slope of the LMC sample. We interpreted this as evidence for 
the existence of an as yet undetected period-dependent system- 
atic error in the IRSB method, but the limited size of the LMC 
sample prevented firm conclusions. 

To put our previous analysis on a firmer basis, we present in 
the present paper new and very accurate radial velocity curves 
for 22 additional LMC Cepheids (Sec|2]i thus almost tripling 
the existing sample for an IRSB analysis to a total of 36 stars. 
To allow a purely differential comparison with solar metallic- 
ity Cepheids we present a similar analysis for 77 Milky Way 
Cepheids in an accompanying paper (Storm et al. 1201 II here- 



Table 1. LMC Cepheids for which we have obtained new radial 
velocity curves with the HARPS and FEROS high-resolution 
spectrographs at ESO-La Silla. 



HV873 
HV876 
HV877 
HV878 
HV881 
HV900 



HV914 
HV1005 
HV1006 
HV1023 
HV2282 
HV2369 



HV2405 
HV2527 
HV2538 
HV2549 
HV5655 
HV6093 



HV12452 
HV12505 
HV12717 
Ul 
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Fig. 1. The photometric V and K light curves, and the (V - K) 
colour curve as used in our IRSB analysis for the LMC Cepheid 
variable HV2282. This is a typical set of photometric light and 
colour curves in the present study based on the optical photome- 
try from OGLE-III (Udalski et al.. 120051 Soszy hski et al. 120081 . 
and the near-IR photometry from Persson et al. ([2004). 

inafter Paper I). In that paper we also combine the new empiri- 
cal constraints on the projection (p-)factor, converting Cepheid 
radial into pulsational velocities, as obtained from the present 
LMC study and the recent HST parallax measurements to ten 
Galactic Cepheids (Benedict et al. 2007), to obtain a new empir- 
ical calibration of the p-factor relation to be used in our IRSB 
distance analysis. 

In Sec[3]We briefly describe the IRSB method in its present 
form and we discuss the choice of the /^-factor relation. We 
present the resulting LMC Period-Luminosity relations in op- 
tical and near-infrared bands in Sec|4]where we also compare to 
the SMC and Milky Way samples to estimate the effect of metal- 
licity on the PL relations. In Sec|5] we discuss our results and 
compare them to other recent investigations. Sec|6] summarizes 
our results. 



2. The Data 

2.1. The LMC Cepheid sample 

We have selected 22 Cepheids in the LMC (see TabQ} show- 
ing a wide range of periods, all having high-quality near-IR 
lightcurves from Persson et al. (2004) and very accurate opti- 
cal photometry from OGLE-III (Udalski et al. 2008 Soszyhski 
et al. I2008I) . From their optical lightcurves, all these Cepheids 
are clearly fundamental mode pulsators. For these stars we have 
obtained radial velocity curves as will be described in the next 
section. A typical data set (HV 2282) is shown in FigJT] 
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Gieren et al. (2005a) performed a first IRSB analysis on a 
sample of 13 Cepheids in the LMC, six of which belonging 
to the young blue, massive cluster NGC1866. These stars will 
all be used in the following analysis. To this sample we have 
added the star HV2827. The radial velocity data for this star 
comes from Imbert et al. ( 1 19871 1, the optical photometry from 
Moffet et al. (1998 ) and the near-IR photometry from Persson 
et al. ( 120041 1. For all the stars we have transformed the near-IR 
photometry to the SAAO system using first the transformations 
given by Persson et al. (2004) to transform the LCO data to the 
CTIO system, and then applying the transformations given by 
Carter (119901 ) to transform from the CTIO to the SAAO system. 
In this way we ensure to be on a common system with the Milky 
Way sample presented in the accompanying Paper I, and to be on 
the same photometric system which was used by Kervella et al. 
(2004 ) for the calibration of the surface-brightness versus colour 
relation. 

2.2. Spectroscopic data 

Using the HARPS instrument on the ESO 3.6m telescope and 
the FEROS instrument on the ESO/MPG 2.2m telescope, both 
at ESO La Silla, Chile, we have obtained 509 new radial veloc- 
ity measurements for the above described sample of 22 LMC 
Cepheids. HARPS was used for 20 of the stars while FEROS 
was used for two additional Cepheids, HV914 and HV12717. 

HARPS (Mayor et al. 2003 ) is a fiber-fed cross-dispersed 
echelle spectrograph mounted inside a vacuum vessel for im- 
proved wavelength stability. It has a spectral resolution of R — 
115000 (3.2 pix) and the wavelength range reaches from 380 
to 690nm. HARPS was used to observe 20 LMC Cepheids in 
the sample. Radial velocity observations were obtained during a 
single, 16-night, visitor observing run in December 2005 which 
were later complemented with service mode observations in the 
following two seasons to cover phase gaps in the velocity curves. 

We used the standard reduction pipeline using cross- 
correlation with a G2-type spectral mask (Baranne et al. 1996 
Pepe et al. 120021 ) to extract radial velocities from the spectral 
data. 

FEROS (Kaufer et al. 1999) is also a fiber-fed crossed- 
dispersed spectrograph. It has a resolution of R - 48000 
(2.2 pix) covering the wavelength range from 360 to 920nm. 
The Cepheid radial velocities were extracted using the cross- 
correlation method with the FEROS pipeline. 

Both instruments are known to be very stable with wave- 
length drifts well below 100 m/s. As our main objective was to 
obtain well covered, accurate radial velocity curves we made rel- 
atively short exposures with resulting low signal to noise ratios 
of 5 to 10. This was found sufficient to keep the errors of the 
individual radial velocity measurements below 100 m/s, which 
meets the precision we desired. 

The Heliocentric Julian dates and the individual radial veloc- 
ities are tabulated in Tab|2]and the radial velocity curves are plot- 
ted in Figj2] A few data points appear to be mis-identifications 
as they are close to the systemic velocity of the star but do not 
fall on the velocity curve (see HV873, HV878, and Ul). They 
have been marked with crosses in the figure and are not consid- 
ered in the further analysis. In the case of HV1005 we found two 
points with similar offsets from the radial velocity curve. These 
two data points were obtained two years later than the majority 
of the data, obtained during the sixteen consecutive nights (see 
above), thus suggesting the possibility of orbital motion for this 
star. We have simply disregarded those two points in the follow- 
ing analysis. 
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Table 2. The new radial velocities for the 22 LMC Cepheids. 
The complete table is available in the electronic form from the 
CDS. 



ID 


HJD 


RV 




(days) 


(kms -1 ) 


HV873 


2453701.7205 


247.46 


HV873 


2453702.7489 


250.14 


HV873 


2453703.7567 


252.64 


HV873 


2453704.7536 


255.15 


HV873 


2453705.7469 


257.52 


HV873 


2453706.7521 


259.96 



3. The IRSB method 

3.1. The surface brightness-colour relation 

The infrared surface brightness (IRSB) method is a variant of the 
Baade-Wesselink method originally developed by Barnes and 
Evans (1 19761 in the optical spectral range. It matches the angular 
diameter variation of a Cepheid as determined from photometry, 
with the radius variation of the pulsating star as determined from 
an integration of its pulsational velocity curve. 

In Paper I we present a more detailed discussion of the 
method, which relates the surface brightness parameter, Fy, to 
a colour index (V - A")o to determine the angular diameter vari- 
ation of the star. 

Here we use the relation 

F v = 4.2207 -O.lVo- 0.5 log 6 (1) 
= -0.1336(V - K) + 3.9530 (2) 

as determined by Kervella et al. (2004). 

3.2. The projection factor 

Knowledge of the projection (p-) factor which converts an ob- 
served radial velocity of a Cepheid into the pulsational velocity 
at its surface is crucial for any kind of Baade-Wesselink analy- 
sis. The ^-factor has to take into account not only the geomet- 
rical projection effects across the observed stellar disk but also 
fold this with the effect of limb darkening and possibly take into 
consideration non-LTE effects due to the dynamic behaviour of 
the pulsating atmosphere. 

Recent empirical (Gieren et al., 12005 al l and theoretical 
(Nardetto et al.. l2007l & l20091 l wo rk has shown that the p-factor 
relation from Hindsley and Bell (1989), which shows a slight 
period dependence of the p-factor and was used in our early 
work (Gieren et al. 1993), is not appropriate. In fact the new 
data presented here (see Sec l4.3b largely supports the empiri- 
cal findings of Gieren et al. (2005a) that the p-factor depends 
quite strongly on pulsational period or the derived distances to 
the LMC Cepheids become period dependent, which is clearly 
unphysical. We have discussed this issue in more detail in Paper 
I. In that paper, we find the best fitting relation to be: 

p = 1.550(+0.04)-0.186(+0.06)logP (3) 

which we will adopt here. 

It is important to note that the exact choice of the /^-factor 
relation has no bearing on the conclusions regarding the effect 
of metallicity on the PL relations in the present paper, as long 
as the ^-factor relation used in the analysis is the same for both 
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Fig. 3. The IRSB-based distance moduli of 36 Cepheids in the 
LMC using the new p-factor relation derived in Paper I, plotted 
as a function of their pulsation period. The moduli have been 
referred to the LMC barycenter using the LMC disk model of 
van der Marel and Cioni (2001). 

samples, the Milky Way and the LMC. Any change in the re- 
lation would affect both Milky Way and LMC Cepheids in the 
same way but the differentials would remain the same. This of 
course is based on the implicit assumption that the /^-factor is in- 
dependent of metallicity. Work is currently underway to further 
investigate this point. 

3.3. The reddening 

The reddenings for the LMC stars were taken from the catalogue 
of Persson et al. ([2004) whenever possible. For the NGC1866 
stars we adopted values of E(B - V) = 0.06 following the dis- 
cussion in Storm et al. (2005). The reddening is very low for 
all the stars and the star-to-star variation is very low as well, so 
reddening errors can only marginally affect the resulting Period- 
Luminosity relations, even in the optical bands. We have adopted 
the reddening law from Cardelli et al. ( 1 19891) with a total-to- 
selective absorption of Ry = 3.23 following the discussion in 
Fouque et al. (2007). In this way we are proceeding exactly in 
the same way as for the Milky Way Cepheids in Paper I. 

3.4. Adopted metallicity 

We adopt the mean metallicity estimates for the Milky Way, 
LMC and SMC Cepheids from the discussion by Luck et al. 
( 11998 ) based on new and literature measures for Cepheids and 
supergiants. They find for the Milky Way [Fe/H] = +0.03 based 
on 69 stars with cr = 0.14, for the LMC [Fe/H] = -0.34 based 
on 32 stars with cr = 0.15 and for the SMC [Fe/H] = -0.68 
based on 25 stars with cr = 0.13. These values are largely 
supported by Romaniello et al. (2008) who found values of 
-0.30 dex and -0.75 dex respectively. We note that the Cepheids 
in NGC1866 might be slightly more metal-poor, [Fe/H] = -0.5, 
based on the spectroscopic study of three cluster giants by Hill 
et al. ( 120001 ) and the photometric measurement for NGC1866 of 
-0.48 + 0. 1 8 by Hilker et al. (119951 from Stromgren photometry, 
but we do not make this distinction in the following. 
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Fig. 2. The new radial velocity curves for 22 LMC Cepheids from HARPS and FEROS (HV 914 and HV 12717) data. The crosses 
indicate points that have not been used in the analysis. The radial velocity range in each panel is 120 km s~' so the amplitudes in 
the panels are directly comparable. The panels are arranged according to increasing pulsational period. 
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Table 5. The derived LMC true distance modulus and the slope 
of the individual Cepheid distance moduli as a function of their 
pulsation period, for different p-factor relations (see text) p = 
P P + a p log P 



"p 


Pp 


(m - M) 


(m — 


M)o slope 


-0.03 


1.39 


18.50 ±0.04 


0.31 


±0.10 


-0.08 


1.31 


18.26 ±0.04 


0.22 


±0.10 


-0.08 


1.455 


18.50 ±0.04 


0.24 


±0.10 


-0.186 


1.550 


18.45 ±0.04 


0.00 


±0.10 



4. Results 

4.1. Distances and absolute magnitudes 

We have applied the IRSB method as described in the previ- 
ous section to find the distances and absolute magnitudes for 
the LMC Cepheids reported in Tab|3] and for the SMC sample 
of Storm et al. d2004l > in Tab|4] The absolute magnitudes are 
intensity-averaged magnitudes. 

4.2. The distances to the LMC and the SMC 

On the basis of the individual distances to the Cepheids in the 
LMC we can now determine the distance to the LMC barycen- 
ter. First we correct the measured distance of each Cepheid 
by the distance offset from the LMC disk model of van der 
Marel & Cioni (2001) as tabulated in the last column of Tab|3] 
and then compute the mean of the resulting values. We find 
(m - M)o(LMC) = 18.45 + 0.04 (random error only) with a stan- 
dard deviation of 0.22 mag. As the Cepheids are well distributed 
across the face of the LMC the uncorrected mean modulus is 
identical to within 0.01 mag to the adopted barycentric distance 
modulus value. 

For the SMC we only have five stars so the random error is 
significantly larger than is the case for the LMC. Additionally, 
the SMC is well known to exhibit more pronounced depth ef- 
fects than the LMC, so some additional scatter is to be expected. 
We find from the five Cepheids an unweighted mean value of 
(m - M)o(SMC) = 18.92 + 0.14 with a standard deviation of 
0.32 mag. The offset in modulus with respect to the LMC of 
0.47 mag is in excellent agreement with the very accurate value 
of 0.44+0.05 mag found by Cioni et al. (2000) from an investiga- 
tion of a large sample of tip of the red-giant branch (TRGB) stars 
in the reddening insensitive near-infrared, thus giving us confi- 
dence that our five stars are indeed representative of the SMC as 
a whole. 

4.3. Constraints on the p -factor relation 

As discussed in Sec l3.2l we can use the IRSB distances to the 
36 LMC Cepheids to place constraints on the p-factor relation. 
We have rederived the distances to the individual LMC Cepheids 
using different p-factor relations, as discussed in more detail in 
Paper I, and we computed the individual LMC barycentric dis- 
tances as before. In Tab|5]we list the resulting mean LMC dis- 
tance moduli as well as the slope of the distance modulus as a 
function of pulsation period, for each of the /^-factor relations 
we tested in this process. The function distance modulus ver- 
sus pulsation period must of course have a zero slope, and with 
the present large sample of 36 stars we are able to place strong 
constrains on the slope. In the table it can be seen that the orig- 
inal Hindsley and Bell (1989) relation gives a very plausible 
LMC distance but a strong and unphysical period dependence 




log(P) 



Fig. 4. The K-band period-luminosity relation for our sample 
of LMC Cepheids based on the absolute magnitudes determined 
with the IRSB method as calibrated in Paper I. The best fitting 
line from the data is overplotted in black, and the best fitting line 
to the Milky Way sample of Paper I is overplotted with a dashed 
line in red. It is appreciated that the LMC and Milky Way PL 
relations agree extremely well both in slope and zero point. 

of the distance moduli. Similarly the Nardetto et al. ( 2009jl re- 
lation gives a slope which deviates at the level of 2cr from the 
uncorrected relation and the LMC distance becomes uncom- 
fortably short at 18.26 mag. The third relation uses the slope of 
the p-factor relation from Nardetto et al. ([2009) but with a zero 
point which yields IRSB distances to Milky Way Cepheids in 
agreement with the recent parallaxes from Benedict et al. (120071 1 . 
The last column corresponds to the relation derived in Paper I 
which we have adopted here, which yields LMC Cepheid dis- 
tance moduli which are totally independent of their pulsation pe- 
riod. This is borne out in Figj3] where we show the LMC center 
distance moduli based on the individual LMC Cepheid moduli 
corrected for the disk model of van der Marel and Cioni (2001 ) 
as a function of log P. It is evident that with the adopted p-factor 
relation there is no significant correlation between pulsational 
period and derived distance to the LMC Cepheids. 

4.4. The LMC Period-Luminosity relations 

In Tab|6]the Period-Luminosity relations we obtain for the LMC 
in the different photometric bands are given in the form: 

M = a(+0-(a))[log(F) - 1 .0] + b(±cr(b)) (4) 

The relations have been determined from a linear regression 
to the absolute magnitudes and log(P) values given in Tab[3] In 
the table the corresponding relations for the Milky Way sample 
determined in Paper I are given for comparison. These relations 
are based on Cepheid distances calculated with exactly the same 
precepts as in this paper, using the same p-factor relation and 
IRSB calibration so we can perform a direct comparison of the 
relations. The relations in the K- and V-bands as well as in the 
V, (V - 1) Wesenheit index are shown in Figs|4]|6]with the Milky 
Way relations overplotted. 

The dispersions of the Milky Way and LMC samples around 
the ridge line PL relations given in Tab|6] are very similar for 
all the bands, and we note that there is not a big difference in 
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Table 3. Distances and absolute magnitudes for the LMC Cepheids calculated using the precepts given in Paper I, using the new 
p-factor relation. In columns 3 and 4 we give the IRSB distance and the associated error estimate. In columns 5 and 6 follows the 
distance modulus, and in the columns 7 through 13, we give the absolute magnitudes in the different bands as well as the Wesenheits 
indices. Finally in column 14 we give the adopted reddening and in column 15 the adopted phase shift between spectroscopic and 
photometric angular diameters. Column 16 gives the magnitude correction to refer the distance moduli to the LMC barycenter. 
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Table 4. Distances and absolute magnitudes for the SMC Cepheids calculated by using the same precepts as in Table 3. Column 
headings are the same as in the previous table. 
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the dispersion as a function of photometric band. The J and K- 
band values are only marginally smaller than those for the op- 
tical bands and the Wyi index performs only slightly better than 
the V-band. If however we consider the LMC Cepheids to be at 
the same distance we can determine apparent PL relations with- 
out involving the individual distances from the IRSB method. In 
this case the dispersion reduces in Wyi from 0.24 to 0.15 mag, 
in V from 0.26 to 0.19 mag and in the ^-band the value even 
decreases from 0.21 to 0.09 mag. Removing only two stars from 
the Wyi relation further reduces that disperion to 0.09 mag il- 
lustrating that the uncertainty on these numbers are quite large 
due to the small sample size. The fact that the dispersion can 
be reduced so much by adopting a common distance to all the 
stars suggests that the observed dispersion in the PL relations is 
dominated by errors in the distance moduli rather than intrinsic 
dispersion in the luminosities or reddening errors and the stan- 
dard error on the individual distance moduli is about 0.2 mag. 



The dispersion of the LMC and MW samples are also very 
similar suggesting that the LMC data quality is equal to that for 
the Milky Way sample. Barnes et al. (2005 ) performed a careful 
Bayesian analysis of the data presented in Storm et al. (2004) and 
found that the formal error estimates from the regression fitting, 
which are also the errors tabulated here in Tabj3] were system- 
atically underestimated. On average these errors should be mul- 
tiplied by a factor of 3.4 to be in agremeent with the Bayesian 
error estimates. In fact the error estimates for the distance mod- 
uli presented in Tabf3] have a mean value of 0.05 1 mag, which 
results in a revised mean error estimate of 0.17 mag, very sim- 
ilar to the standard error of 0.2 mag estimated above from the 
^-band PL relation. 

4.5. The effect of metallicity 

The LMC Cepheid PL relations in Tab|6] can be compared di- 
rectly to the relations derived in Paper I for the Milky Way 
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Fig. 8. The zero point offsets, AM = M - M com b, for each band with respect to the reference PL relations from Paper I as a function 
of sample metallicity. The full lines show the weighted linear regression fit to the data. It can be seen that especially in the limited 
but for extragalactic distance determination most important metallicity range from LMC to solar abundance, the zero point offsets 
are in general indistinguishable from zero. 
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Fig. 5. The V-band PL relation for our sample of LMC Cepheids 
based on the absolute magnitudes determined by the IRSB 
method. Overplotted is the linear regression fit (solid line) as 
well as the corresponding Milky Way relation (dashed line). 

sample which have been tabulated for convenience in Tabj6] 
The slopes of the relations are in excellent agreement in the 
case of the near-infrared J and K bands. In the V band and in 
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Fig. 6. The optical (V - I) Wesenheit index PL relation for 
our sample of LMC Cepheids, based on the absolute magni- 
tudes determined by the IRSB method as calibrated in Paper I. 
Overplotted is the linear regression fit (solid line) as well as the 
corresponding Milky Way relation (dashed line). 

the Wesenheit indices, the agreement is slightly worse but well 
within one <x, whereas the / band slopes differ at the level of 
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Fig. 7. The /T-band absolute magnitudes for our sample of SMC 
Cepheids (filled triangles) overplotted on the LMC Cepheid PL 
relation (FiglUi, crosses are the individual LMC Cepheids. 

Table 6. The Period-Luminosity relations for the LMC in the 
various optical and near-infrared bands as determined from a lin- 
ear regression to the absolute magnitudes from the IRSB analy- 
sis. The relations have been fitted in the form M - a[log(f) - 
1] +b. The fourth column gives the dispersion around the fit. Our 
results for the Milky Way Cepheids are included for comparison. 
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Table 7. The offsets of each metallicity sample (Milky Way, 
LMC and SMC) with respect to the reference PL relations from 
Paper I, i.e. assuming a fixed PL relation slope irrespective of 
metallicity. The last column gives the slope, y, of the weighted 
linear regression fit to the data. 
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about one cr. Since the slopes of the PL relations for both LMC 
and Milky Way Cepheids are very well constrained from our 
samples, we conclude that our results present strong evidence 
that the slope of the Cepheid PL relation, particularly in the 
near-infrared J and K bands, is identical for the solar metallicity 
Milky Way and more metal-poor LMC samples. These PL rela- 
tions are thus independent of metallicity in the range from solar 
to LMC metallicity. The universality of the PL relation slopes 
appears to be confirmed in this metallicity range. 

We can further extend the metallicity baseline by including 
the 5 SMC Cepheids in the sample. Given the low number of 
Cepheids and the narrow range of periods for these stars, we 
cannot constrain the slope of the PL relation at this metallicty, 
but we can constrain the zero point offset. In FigQwe have over- 
plotted the SMC Cepheids on the LMC Cepheids for the K-band 
and obviously there is good agreement. Considering the excel- 
lent agreement between the LMC and the Milky Way from Tab(6] 
this further supports the universality of the zero point of the K- 
band relation. 

We can quantify any offsets in the PL relation magnitude 
zero points as a function of metallicity by comparing the zero 
point offset for each of the three samples with the reference 
PL relations determined in Paper I on the basis of all available 
Cepheids. For each band and sample we computed the mean 
magnitude offsets, tabulated them in Tab|7] and plotted them as 
a function of metallicity in FigJS] For each band we have fit the 
weighted least square regression line to the magnitude offset as 
a function of metallicity. The resulting metallicity effect slopes 
have been tabulated in Tab(7] We estimate the error on the slopes 
to be of the order 0. 10 mag dex 'based on the error bars on the 
individual points and the fact that the baseline is smaller than one 
dex. From the FigJS] and TabQit can be seen that a zero metal- 
licity effect cannot be ruled out by the data in most cases but 
that the effect is significant and the strongest for the Wvi index. 
For extra-galactic distance determination most Cepheid samples 
will have metallicities in the range from LMC to solar and the 
zero-point offsets in this range is clearly very small in all the 
bands. 

The only bands for which we have excellent agreement be- 
tween the PL relation slopes are the J and K bands, and the zero- 
point offsets in both of these bands are small so they each form 
an excellent basis for a standard candle. In the V band and the 
Wjk index the PL relation slopes are marginally in agreement 
and also for these two bands the zero-point offsets are small. 
The / band and the Wvi index have significant differences in the 
PL slopes. The Wvi index zero-point is also significantly more 
metallicity dependent than is the case for the other bands mak- 
ing this index inferior to the J and /T-bands as a standard candle. 

We note that we have assumed a simple linear metallicity 
dependence, but it is of course entirely possible that the depen- 
dence, if present, has a more complex functional form. 

We emphasize that the metallicity effect on both the slopes 
and zero points discussed here is entirely independent of the 
choice of the p-factor relation, whereas the absolute values of 
the slopes and zero points as well as the resulting LMC and SMC 
distances do depend on the adopted /^-factor relation. 
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5. Discussion 

We can now compare the slopes of our LMC PL relations with 
independent measures from apparant magnitudes versus log(P). 
In the near-infrared we use the relations from Persson et al. 
d20041 > which is based on 92 LMC Cepheids with excellent / 
and K light curves. They found slopes of -3.261 + 0.042 and 
-3.153 + 0.051 in K and J respectively, in excellent agreement 
with our values of -3.28 + 0.09 and -3.22 + 0.09. In the optical 
we can compare to the relations from the OGLE project (Udalski 
2000) and we similarly find very good agreement. They find in 
the V and /-bands slopes of -2.775 + 0.031 and -2.977 + 0.021, 
to be compared with our values of -2.78 + 0.1 1 and -3.02 + 0.10 
respectively. The fact that we reproduce the slopes of the PL re- 
lations in the LMC so well strongly supports our empirical cali- 
bration of the p-factor used with the IRSB method and confirms 
our earlier results (Gieren et al. 2005a) based on a much smaller 
sample of LMC Cepheids. 

The LMC distance modulus of 18.45 + 0.04 is only slightly 
shorter than the "canonical" value of 18.50 + 0.10 as adopted 
by Freedman et al. (120011 ) and by the ARAUCARIA project 
(Gieren et al. 2005b, Pietrzyriski et al. 2010a), and slightly 
longer than the value of 18.39 ±0.01(random) ±0.07(systematic) 
suggested by Freedman and Madore (2010) in their recent re- 
view. It is thus in good agreement with most recent results, as 
would be expected considering that we (see Paper I) have cal- 
ibrated the IRSB method to match the distances to nine Milky 
Way Cepheids with direct parallaxes from Benedict et al. (|2007 ). 
These distances are accurate to 4-10%, which seems the largest 
single remaining systematic uncertainty in the IRSB method as- 
suming that the p-factor relation has been accurately established 
from our new and improved constraints (which is supported by 
the arguments given above). The total systematic uncertainty of 
our present determination of the LMC distance is difficult to 
quantify, but expected to be about ±5%, which is better than 
most other techniques which have been used to measure the 
LMC distance over the years. An exciting possibility to con- 
strain the IRSB method even further, in a very direct way, will 
be the comparison of the IRSB-determined distance to the LMC 
Cepheid OGLE-LMC-CEP0227 which is a member in a double- 
lined eclipsing binary system (Pietrzynski et al. 2010b). Work is 
underway to determine the distance to this binary from orbital 
analysis. 

In Sec l4.5l we compared the slopes of the PL relations from 
the Milky Way and LMC samples finding insignificant differ- 
ences in the near-IR bands and small (but possibly also non- 
significant) differences in the optical bands. This result is in good 
agreement with the findings by Gieren et al. (2005a) but at odds 
with earlier findings by us (Fouque, Storm and Gieren, 2003 
Storm et al. 120041 ) where we found that the slope of the galactic 
relations were significantly steeper than those in the LMC and 
SMC. The reason for this difference was the inappropriate p- 
factor relation used in the earlier work for determining the Milky 
Way PL relations. Note that in the present work the choice of 
the p-factor relation has no bearing on the difference in slope 
between the Milky Way and LMC PL relations as we are now 
applying the method to both samples of Cepheids, this was not 
possible previously as the necessary data for the LMC stars was 
not yet available. Sandage, Tammann, and Reindel (2004) also 
argued that there is a significant difference in slope between the 
Milky Way and LMC PL relations in the optical bands, the Milky 
Way relation being steeper. This conclusion was largely based on 
these older IRSB results, as well as on Cepheids in OB associa- 
tions from Feast ( 119991 ). As we show in Paper I the revised IRSB 



distances are, apart from a small zero point offset, in very good 
agreement with the latest results on open cluster Cepheids from 
Turner (2010) and does not exhibit a period dependence. 

We do not support the conclusion of a strong effect of metal- 
licity on the slope of optical (V, /) Cepheid PL relations as 
reached, for example, by Tammann et al. d201 II ) or Storm et al. 
(120041 ). in fact we find that for these bands the Milky Way PL 
relations are slightly shallower than the LMC relations, only for 
the Wjk index we find that the Milky Way relation is steeper than 
the LMC relation. 

In Tab|7] we summarized the PL relation zero point varia- 
tion as a function of metallicity. We have made the simplest as- 
sumption of a linear relation and found the slopes of the zero 
point offset versus metallicity relations, y. Clearly the effects 
in the near-IR are small and even in the V band the effect ap- 
pears to be small, albeit with the opposite sign as for the other 
bands. The largest effect we find is in the Wyi index and our 
value of y(Wvi ) = -0.23 + 0.10 mag dex~Hs in excellent agree- 
ment with the value of -0.24 ± 0.05 mag dex'found by Sakai et 
al. (|2p04) and other measurements which have been adopted in 
the recent review paper by Freedman and Madore (2010). This 
result is slightly at odds with the findings of Bono et al. (2010) 
who on the basis of data on 48 external galaxies finds no sig- 
nificant metallicity effect on neither slope nor zero point for the 
Wesenheits indices Wvi and Wjk ■ At the same time we do agree 
with them that the PL relation slopes are less affected in the J 
and -bands and more affected in the optical V and / bands and 
we also agree on the most likely sign of the effect namely that 
metal-rich PL relations are shallower than metal-poor ones. 

It is interesting to note that the most significant metallic- 
ity effect is found for the Wy/ index and that this effect can 
largely be attributed to the colour difference between Milky Way 
and LMC Cepheids as found by Sandage, Tamman and Reindl 
(2004). They compared period-colour relations for both Milky 
Way and LMC Cepheids and found that the LMC Cepheids for 
a given period are bluer by about 0.05 mag. Feeding this into the 
Wesenheits index (Wyi = V — 2.54(V - /)), we find an offset 
of 0.00 - 2.54 x -0.05 = 0.13 mag in the case where the offset 
is fully in the /-band and -0.05 - 2.54 x -0.05 = 0.08 mag in 
the case where the effect is fully in the V-band. Disregarding 
the subtle effects of differences in the slopes of the period- 
colour relations, these results are comparable to the offset of 
0.09 mag, metal-poor Cepheids being fainter, which we find be- 
tween the Milky Way and LMC Wyi relations in Tab|6]at a period 
of 10 days. 

The emerging conclusion based on our data and analysis is 
that for accurate distance measurements to galaxies the /^-band 
Cepheid PL relation is the best suited tool: it is metallicity- 
independent both regarding the slope and the zero point, it is 
very insensitive to reddening, and it has a smaller intrinsic dis- 
persion than any optical PL relation. It is likely, as indicated in 
recent work from Spitzer data, that mid-infrared Cepheid PL re- 
lations are even superior to their near-infrared relations because 
of their even lower sensitivity to reddening, and lower intrin- 
sic dispersions (Madore et al. 120091 ). Yet, their dependence on 
metallicity has still to be investigated and they cannot be ex- 
ploited from the ground making them exceedingly expensive to 
use. 



6. Conclusions and Summary 

We have obtained new and very accurate radial velocity curves 
for 22 LMC Cepheids thus expanding the sample of LMC 
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Cepheids to 36 for which we can apply the IRSB distance anal- 
ysis. We have applied the newly calibrated IRSB technique of 
Paper I to these 36 LMC Cepheids as well as to 5 SMC Cepheids. 
The IRSB analysis yields individual distances from which we 
calculate absolute magnitudes in optical and near-infrared bands. 
These magnitudes define tight Period-Luminosity relations in 
the V, I, J, K bands and in the Wesenheit indices. 

We show that the PL relations are in excellent agreement 
with the observed apparent magnitude versus log(P) relations in 
both the optical and near-IR bands lending strong support to the 
empirical calibration of the p-factor relation in Paper I. 

By comparing the LMC Cepheid PL relations to their Milky 
Way counterparts reported in Paper I that were established with 
exactly the same precepts, we find practically identical Milky 
Way and LMC PL relation slopes in the near-infrared bands, 
arguing for the universality of the Cepheid PL relation in this 
spectral range. The zero points exhibit a slight metallicity effect, 
j(Mk) - -0.10 ± 0.10 mag dex'in the sense that metal-poor 
Cepheids are fainter than metal-rich Cepheids. If we restrict our- 
selves only to the metallicity range between solar and LMC, our 
results are consistent with universal PL relations (in both slope 
and zero point) in the near-infrared J and K bands. 

In the optical bands, we find that the slopes depend weakly 
on metallicity, the Milky Way slopes being slightly shallower 
than the LMC slopes, but this difference might not be significant. 

The optical PL relation zero points exhibit a metallicity effect 
of similar size as in the near-IR albeit the sign is opposite in the 
V-band. The Wyi index shows the strongest zero point effect on 
metallicity where we find y{Wvi ) = -0.23 + 0.10 magdex -1 . 
We stress that the zero-point offsets are based on a rather long 
metallicity baseline ranging from SMC to solar metallicity. 

Our direct distances to the LMC Cepheids leads to a true 
LMC distance modulus of 18.45 + 0.04 mag, with an estimated 
systematic uncertainty of 5% which mainly comes from the un- 
certainty on the HST parallaxes of nine Milky Way Cepheids 
that have been used to define the absolute zero point of the IRSB 
technique. Our distances to the SMC Cepheids leads to an SMC 
distance of 18.92 + 0.14. 

Both the Wyi metallicity effect and the LMC distance which 
we find are in agreement with the latest values adopted by 
Freedman and Madore (120 1 Oi l in their recent review on the 
Hubble constant. 

Considering the significant metallicity effect on the 
Wvi index, not only on the zero-point but most likely also on 
the slope, we argue that the best standard candle is presently 
provided by the PL relation in the /T-band as it is metallicity in- 
sensitive, reddening insensitive and exhibits the lowest intrinsic 
scatter. 
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